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ABSTRACT 

Continuum observations at 350 fim are presented of seven nearby elliptical 
galaxies, for which CO-gas disks have recently been resolved with interferometry 
mapping. These SHARCII mapping results provide the first clearly resolved 
far-infrared(FIR) to submillimeter(submm) continuum emission from cold dust 
(with temperatures 32 K > T > 22 K) of any elliptical galaxy at a distance 
> 40 Mpc. The measured FIR excess shows that the most likely and dominant 
heating source of this dust is not dilute stellar radiation or cooling flows, but 
rather star-formation, that could have been triggered by an accretion or merger 
event and fueled by dust-rich material that has settled in a dense region co- 
spatial with the central CO-gas disks. The dust is detected even in two cluster 
ellipticals that are deficient in HI, showing that, unlike the HI, cold dust and CO 
in ellipticals can survive in the presence of hot X-ray gas, even in galaxy clusters. 
No dust cooler than 20 K, either distributed outside the CO disks, or co-spatial 
with and heated by the entire dilute stellar optical galaxy (or very extended H I), 
is currently evident. 
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Introduction 



Cold gas and dust in nearby elliptical galaxies were detected on ly fairly recently (jSadler fc Gerhard 



19851 : iKnapp et al.lll989l : iLees et al. 



1991c IWiklind fc Henkellll995l ). Compared to that in spi- 



ral galaxies, the cold Interstellar Medium (ISM) in ellipticals is present in relatively small 
amounts and is seen in only 50% to 80% of nearby ellipticals. The source and content of the 
cold ISM in these galaxies is still uncer tain, with optical and far-in frared (FIR) dust-mass 
estimates differing by ~ 10 to 100 (e.g.. iGoudfrooij fc de Jonglll995l ). 



There is some evidence for an internal origin of the dust arising in circumstellar envelopes 
of red giant stars. This internal origin is supported by the detection of ~ 10 /im emission 
in excess of the stellar emission, with the same de Vaucoul eurs profile as the optical stellar 



images and near- to mid-IR photospheric emission (e.g., lAthey et al.l 120021 ) . Some have 



suggested that this dust would come in contact with the hot ISM (~ 10 7 K) and be erode d 



in therm al collisions with ions and dest royed in only ~ 10 7 to 10 8 years (e.g.. lSchweizerlll998l ). 



However, iMathews fc Brighentil (120031 ) have argued that the dust may induce efficient cooling 
of the hot gas on a timescale shorter than the destruction lifetimes of the dust in the hot gas, 
allowing some dust to survive and settle primarily in the centers of the galaxies while the hot 
gas cools. They propose that if this dust-induced cooling is the dominant evolutionary path 
of the dust from stars, the thermal energy released by the hot ISM should be re-radiated 
by dust in the FIR and could be responsible for the central dusty disks and clouds seen in 
ellipticals. 

It has also been suggested that mergers or accretion events that include a gas-rich 
galaxy may be th e dominant source of dust in ellipticals. The total far-IR l uminosities (e.g., 
Temi et al.ll2004j ) and cold gas content (IKnapp et al.lll985c ILees et al.lll99ll ) of ellipticals are 
observed to be uncorrelated with their optical luminosities. This lack of correlation, together 
with HST detections of irregular dust lanes and filamen ts with random orientations with 
respect to the optical major axes (e.g., iTran et al.l 120011 ) is usually interpreted as evidence 
that the dust has been externally acquired. The external origin is likely to be accretion via 
the interaction of an elliptical galaxy with a gas-rich spiral, or in mergers of gas-rich spirals 
that evolve through an FIR-bright (starb urst) phase as they relax to form ellipticals (e.g., 
Toomre fc Toomre 1972 ; Schweizer 19981 ). In fact, if elliptical galaxies are formed through 
spiral mergers, then starbursts are necess ary to produce the high p hase-space densities seen 
in the centers of elliptical galaxie s (e.g. , iRothberg fc Joseph! 12004 . and references therein). 
Numerical simulations by iBarned (120021 ) have shown that 20% to 60% of the gas from two 
spirals merging can form relaxed disks of sizes up to 20 kpc, with some of the gas maybe 
losing angular momentum, falling to the nucleus of the galaxy, and perhaps leading to star- 
formation. 
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Observations of dust emission from elliptical galaxie s have been made using IRAS or 



ISO primarily (e.g.. lKnapp et al.lll989l ; iTemi et al.ll2004l ). These measurements favored the 



warm dust components of the ISM and produced unresolved images of the dust emission. The 
IRAS and ISO derived SEDs of these galaxies clearly show that there are multi-temperature 
components to the dust. Indeed, the presence of an additional, low-level, very cold (< 20 K) 
dust component in ellipticals has been suggest ed by limited detections of unresolved mm 



and submm emiss i on in some nearb y galaxies (IKnapp &: Patten! Il99ll ; IWiklind et al.l 11995 



Leeuw et al.ll2004 ITemi et al.l 12004 ). as well as clearly extended FIR to submm detections 
from the dust lane in the nearest giant elliptical galaxy NGC 5128 that hosts t he active galac- 



tic nucleus (AGN) CentaurusA (e.g.. iLeeuw et al.l 120021 ; iQuillen et al.ll2006l ). Observations 



of CO emission from a nu mber of elliptical s show that many of the m contain significant 
amounts of molecular gas (ILees et al.l Il99ll ; IWiklind fc Henkell 1 19951 ). but many of these 
observations do not spatially resolve the CO emission. 



Recently, interferometric measurements by I Young (J2002J, |2005|) and lWiklind et al.l (119971 ) 
have spatially resolved the CO (1-0) emission in seven elliptical galaxies with spatial resolu- 
tions of 7 and 3 arcsec, respectively. These data imply that the molecular gas is contained 
mostly in rotating disks lying within the central 5kpc of the galaxies. This prompted the 
study by our group to spatially resolve the dust emission in these seven elliptical galaxies. 
Our goals were (a) to determine the physical properties of the dust coincident with the 
molecular gas; (b) to see if there are dust associations with other spatial distributions (e.g., 
stellar distribution); and (c) to determine the most likely heating source (or sources) for 
the dust. We wanted also to address the question of whether these observed gas disks are 
the result of external merger/accretion events, as suggested by the optical morphological 
signatures of the galaxies, or internal mass-losing stars that constitute their elliptical optical 
components. 

The seven galaxies in our sample (NGC 83, N GC 759, NGC807, UGC15 03, NGC 3656, 
NGC4476, and NGC 5666) were all classified bv IWiklind fc Henkell Jl995h as ellipticals, 
based on the determination that their luminosity profiles follow the i? -1 / 4 de Vaucouleurs 
law. However, morphologically the sample includes galaxies from different extragalactic 
environments and represents a spread of possible merger /accretion tracers or ages, from 
galaxies that have been classified as on-going or early-age major mergers (e. g., NGC 365 6, 



Balcells et al.ll2001l) to very-late accretion or quiescent systems (e.g., NG C 807. iMurrav et al. 



2000). Recent CCD imaging of some of the galaxies, e.g., NGC 5666 by lDonzelli fc Davoust 



(120031 ) . have disputed the classification as elli pticals in fav or of SO galaxies, while new 
optical observations of others, e.g., NGC 83 by lYound (120051 ). have confirmed the elliptical 
classification. 
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The sample's updated environments, general properties, and molecular gas content are 
summarized in Table 1. The observational merger or accretion signatures of the sample 
are worthy of note as they fortuitously allow an evolutionary analysis in the context of 
merger /accretion stages and cold dust and gas properties. These merger /accretion signatures 
are described further with other observational characteristics of the galaxies in Section [3] 
below. 



2. Mapping Observations using SHARCII at 350/um 

We observed NGC83, NGC807, UGC 1503, and NGC759 with SHARCII on October 
11-14, 2005, and NGC3656, NGC4476, and NGC5666 on April 23 and 29, 2006. SHARCII 
has a 12 x 32 bolometer array with an instantaneous field of view of 148" x 56". We achieved 
background subtraction and mapping coverage for the galaxies by having the telescope ex- 
ecute a Lissajous scanning pattern with amplitude 60" in azimuth and 40" in elevation. 
During the observations, we maintained a log of air mass and the zenith water vapor; the 
telescope dish-segments were kept in optical alignment using the CSO Dish Surface Opti- 
mization System (DSOS); and pointing calibrations were executed every hour, allowing us 
to achieve a FWHM beamsize of 8.5" and pointing accuracy of about 2" rms. The zenith at- 
mospheric opacity at 225 GHz, during the observations and as measured by the CSO facility 
radiometer, varied from 0.04 to 0.06. 

For flux calibration, we observed either Callisto, Uranus, CRL618, CRL2688, Vesta, 
Arp 220, CIT6, Neptune, or o Cet hourly. The assumed fluxes for the exosolar calibrators are 
10.2 Jy/beam for Arp 220, 19.4 Jy/beam for CRL618, 2.33 Jy /beam for o Cet, 41.6 Jy/beam 
for CRL2688, and 2.42 Jy/beam for CIT6, and, for the solar-system calibrators, in October 
232 Jy/beam for Uranus, 90.0 Jy/beam for Neptune, 5.77 Jy/beam for Vesta, and in April 
95.8 Jy/beam for Callisto. The absolute calibration uncertainty is estimated to be 20%. 

The SHARC II contour maps are shown in Figures [T] to [7] and are the result of about 
3.0 hours of observing on NGC 83, 7.4 hours on NGC 807, 6.3 hours on UCG 1503, 2.3 hours on 
NGC759, 1.9 hours on NGC 3656, 2.4 hours on NCG4476, and 1.3 hours on NGC 5666. The 
data were analyzed using the SHARC II reduction software SHARCSOLVE, which iteratively 
solves for the source image, atmospheric background, and detector gains and offsets. During 
the iterations, a constraint was applied in the map domain: outside a diameter of 37" or 
46", the source image was assumed to be zero. The 46" diameter was applied to the maps 
of NGC 807 and UGC 1503, and the 37" diameter applied to those of the other sources. The 
positional accuracies of the galaxies are based on the pointing observations of the calibration 
sources. However, to optimize the effective resolution, the coadded images are composed 
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of hour-long sub-images whose centroids have been aligned. This relative alignment was 
performed for all objects except NGC 807 and UGC 1503, which were too faint and extended 
to have well-defined centroids after one-hour integrations. The 350 fim fluxes of these maps 
are given in Table [2J The integrated fluxes were obtained in a 46" aperture for NGC 807 and 
UGC 1503 and a 37" aperture for all the other target objects. 



3. Submm Morphology versus Observations from other Wavebands 



The goal of this submm wavelength study is to directly map the emission from cold dust 
that may be c o -spat ial w ith the optica l dust features and/or the CO gas-disks presented by 
Wiklind et al.l (119971 ) and I Young] (120021 . 120051 ) for the sample galaxies, as well as any emission 
from dust not detected by the optical and CO studies. The presented submm mapping is 
sensitive to continuum emission from optically-thin dust and, thus, the observations will 
detect emission from dust that may be deep or in the far-side of the galaxies, and thus not 
easily visible in the optical images. The detected submm distribution and fluxes are analysed 
to assess their galactic associations and the possible mechanisms heating the dust; therefore, 
the submm mapping spatially and physically probes the dust's origin, in complement with 
observations at other wavelengths. 

Below, the distributions of the mapped submm emission are compared to those of galac- 
tic components observed at other wavelengths. Later in Section [5l the observed submm fluxes 
are used to calculate and analyse the heating of dust by star-formation and other possible 
submm-emitting, galactic sources. 

Figures [TJ to [7J show contours of the detected 35 (im dust con t inuum for t he sample and 



the to tal integrated CO (1-0) intensity maps from lYoungi (120021 . 120051 ) and IWiklind et al. 



(119971 ). overlayed on the optical Sloan Digital Sky Survey (SDSS) maps, and where there are 
no SDSS images available on the Digitized Sky Survey 2 (DSS2) images. The SDSS images 
have more reliable photometry than the DSS2 images; and, the SDSS images are plotted 
convolved to a 10" beam similar to that of the CSO dust emission images, while the DSS2 
images are plotted at their full resolution. Figures [fl [2j and [7J also show SDSS i- minus 
g-band color or extinction maps and their linear contours in an attempt to show some of the 
dust's geometry. For a smooth presentation, these color maps are shown convolved to a 2" 
beam. 

The SDSS and DSS2 intensity contours are percentages of the respective peaks in the 
maps and are intended to display areas that may be co-spatial at the different wavelengths. 
The 50%-contour-level diameters determined by eye from these plots after deconvolution are 
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listed in Table [3) For NGC3656, the optical 50%-contour-level diameter is well off the maps 
ffotted in Figure (H and therefore only the optical 70%-contour-level diameter is listed in 
Table [21 These tabulated image sizes were derived after deconvolving beamsizes of 7" for 
the CO images (except for the size of NGC759, where a 3" beam was used) and 10" for the 
submm and SDSS images. For the submm and SDSS images, the diameter measurements 
done by eye were compared with FWHM diameters from 2D gaussian fits assuming elliptical 
shapes and agreed to within about ~ 5%. For the DSS2 images, the diameters listed are 
from the full resolution maps. The relative percentage contours of the optical maps will be 
affected by the distribution of dust in these galaxies, therefore the optical intensity contours 
on those maps are not a simple reflection of the stellar distributions as shown on those maps 
and listed in Table (3) but rather, they show the stars plus their attentuation by internal 
galactic dust of varying distributions, and for the DSS2 images may also show saturation 
from the original photographic images. 

The 350 //m results show that the detected dust continuum of all the presented ellipti- 
cals generally follows the extended CO emission that has been previously detected in these 
galaxies. This is therefore consistent with emission from cold dust that occurs in dense, 
gas-rich star-formation regions, as is comm only seen in mapping ob servations of the Milky 



Way and other star-forming galaxies (e.g., ISeaquist fc Clarkl l200ll ) and confirms CO-FIR 



correlations of unresolved data f rom elliptical gala xies that have suggested that cold gas and 



dust co-exist in ellipticals (e.g.. iLees et al.lll99ll ). Individual galaxies are discussed in the 
following sub-sections. 



3.1. The Early-Age Major-Merger/ Accretion Elliptical NGC 3656 



Resolved imaging of NGC 3656 by lBalcells et al.l (120011 ) showed that H I occurs in shells 
and tidal tails that are thought to be signatures of a recent (or early-age), major merger of 
this galaxy. The HI is more extended than the galaxy's main optical c omponen t . Sen sitive 
CCD optical images of th e dust features in NGC 3656 were studied by iBalcellsl ( 119971 ) and 
Rothberg fc Joseph! (120061 ). who suggest that the dust is a remnant of a dust-rich merger 
that involved this galaxy's primary optical elliptical component. 

Figured] shows that large-scale, 350 /im-dust continuum spurs in the south east of NGC 36 56 



and 


Balcellsl f 


1997) 


Balcellsl ( 


1997 


) in t 



The 350-/zm continuum spurs also follow the optical spurs shown by 



tures of a warped-disk of dust that has been reported in the center of NGC 3656 (e.g. 
19971 ) and described further below. 



Balcells 
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Figure [T] also shows that the extinction contours of NGC 3656 are very asymmetric, 
depicting strong extinction in the galactic east and dense north-south central dust distribu- 
tion. The dust in NGC 3656 is probably in a thin, warped disk that provides the strongest 
extinction where the disk is tangent to the line of slight. In the center, the thin disk is there- 
fore tangent to the line of sight and thus provides the strong optical extinction and submm 
emission seen there, while it is more face-on or in the background in the western part of the 
galaxy. In the east of the galaxy, the disk is either more edge-on or in the foreground, and 
thus exhibits the moderate optical extinction seen there (see Figured]). 

A thin, warped disk similar to that d escribed for NGC 3656 has also been o bserved in 
another early- merger galaxy CentaurusA (ILeeuw et al.ll2002l ; iQuillen et al.l 120061 ) . showing 
similar extinction and submm emission features as seen in NGC 3656. In CentaurusA, it has 
been speculated that the disk has formed from dust-rich ma terial accreted from a merger 
that involved a gas- rich spiral galaxy (e .g. . IQuillen et al.ll200q) . The similarity of the warped 
dust-disk together with the HI shells (IBalcells et al.l 120011 ) detected in NGC 3656 to those 
same features observed in Centaurus A suggests these two galaxies have a similar evolutionary 
history and that NGC 3656 may indeed be a younger, more distant analog of CentaurusA. 



3.2. The Intermediate-Age Merger/ Accretion Elliptical NGC 5666 



The HI in NGC 5666 extends well beyond th e optical galaxy; and, it is in circular 
orbits that indicate it is in dynamic equilibirum (ILake et al.l 119871 ). since a not too re- 
cent (or int ermediate-age) merger /ac cretion event that is thought to have occured in this 
galaxy (e.g., 



Donzel 



i fc Davoust 



20031 ) . Using the Canada France Hawaii Telescope (CFHT), 



Donzelli fc Davoustl (120031 ) obtained optical CCD spectroscopy and imaging that revealed a 
"spiral" dusty structure with tidal tails and H II regions this galaxy. The "spiral" and tidal 
tails are revealed after subtracting a model of the bulge plus disk of the host galaxy. 

Figure [2] shows that CO emission and submm continuum contours of NGC 5666 follow 
each other in the high-brightness region greater than the 50% contour level. Outside the 
10"-radius the CO emission contours are more North-South while the submm ones (which 
are not detected at a high signal-to-noise) are East- West and slightly more extended than 
the CO. In the eastern part of the galaxy, the submm and optical SDSS i-band contour levels 
follow each other down to the 30% level. However, in the western part, these contour only 
follow each other to the 70% levels. The asymmetry in the contour correspondence supports 
the idea that the dust distribution or extinction in NGC 5666 is not even, but for example 
in a "spiral" structure as suggested by iDonzelli fc Davoustl (120031 ) . 
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Figure [2] also shows a red extinction region within the 10" diameter that coincides with 
the CO emission and sub mm continuum peaks and the location of the "spiral" that was 
revealed in this galaxy by iDonzelli &: Davoustl (120031 ). This extinction map also shows a 
patchy region outside 15"-radius that may coincide with foreground dust patches or the 
location of the optical tidal tails connected to the "spiral" dusty structure. 



3.3. The "Late-Merger/ Accretion" Ellipticals NGC 83 and NGC 759 



NGC083 and NG C 759 are elliptical galaxies with surface brightness profiles that are 
fit by the r 1 / 4 law (e.g. lYoungll2005l ; IWiklind et al.lll997l ). Using deep CCD V and R images, 
Young ( 20051 ) confirmed in detail that NGC 83 is classically fit by the r 1 / 4 law from about 
radii 6" to 80" and that within about radii 6" it has a regular-shap ed cold gaseous dust disk. 
NGC 759 is also an elliptical galaxy in which IWiklind et al.l (119971 ) earlier showed that there 
is a compact, star-forming, central cold gaseous disk 

Consistent with dust emitting c omponents t hat are likely co -incident with the co mpact 
gaseous disk detections in NGC 83 (TYoungl l2005 h and NGC 759 Jwiklind et alJll997h . both 
these galaxies are not resolved in the submm images presented in Figures [3] and 0J The 
DSS2 blue and re d optical imag es of these galaxies are saturated; and here, the V — R 
extinction maps of lYoungl (120051 ) are compared with SHARCII emission maps of the dust. 
The compact size derived from the V — R extinction map of the dust disk in NGC 83 by 
Young] (120051 ) is smaller in extent than the 350 //m continuum beam, consistent with the 
deconvolved, compact submm size of this galaxy measured from the presented data, and 
that submm emission from the dust-disk of NGC 83 (or NGC 759) being simply unresolvable 
with SHARC II. 



The central disks in NGC 083 and NGC 759 could represent a late stage in the merging 
of two disk galaxies or accretion of gas-rich material. The gas could have lost momentum in 
the merging process (or since the accretion) and fallen to the center of the galaxy. If this gas 
component is capable of forming stars, it could produce the high phase density characteristic 
of ellipticals. 



3.4. The Quiescent or Very-Late- Accretion Ellipticals UGC 1503 and NGC 807 

Of the ellipticals presented here, UGC 1503 and NGC 807 are the only object that are 
(1) field galaxies and (2) have no currently identified optical morphological signatures of a 
merger or accretion event. Therefore, they are referred to as quiescent or very-late-accretion, 
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to mean that if there has been any possible accretion event during their evolutionary lives, 
that happened a very long time ago, and the galaxies have now settled beyond any such event. 
Being field galaxies also means that there is no near neighbor that may have interacted with 
them in the near past. 

UGC 1503 and NGC 807 have intrinsic 50% contours of the CO-emission and 350 /im- 
continuum fluxes that are comparable in extent to and only slightly smaller than those of 
the optical. These galaxies probably have dust that is distributed in slightly edge-on disks 
with their ne ar-side in the northeast and northwest for UGC 1503 and NGC 807 respectively. 



Young! (120021 ) report that in NGC 807, 70% of the CO emission is detected in the south of 



the galaxy (see Figures [6]), suggesting the bulk of the disk's content is in the southern parts. 



3.5. The Accretion Virgo-Cluster Dwarf Elliptical NGC 4476 



Despite dir ected searches, H I 



NGC 759) (e.g., 



flLees et al.lll991 



ras not been detected in the cluster galaxy NGC 4476 (or 



Lucero et al.l 120051); this i s in spite of CO-gas detections in these galaxies 



Wiklind fc Henkellll995l : IWiklind et al.l 1 199 71 : lYound 120021 ) . It is thought 



the HI may been destroyed by the h ot cluster gas or, in the case or NGC 4476, by ram- 
pressure stripping (ILucero et al.l 120051 ) . 



Using HST V (F555W) imaging. Ivan Dokkum & Franxl (119951 ) showed that NGC 4476 
has a regular, "spiral-like" dust lane that is about 20" in diameter. They determined a 
difference of 19 degrees between the major axis position angle of the dust lane and that of 
the stellar component and used the evidence in this and a sample or other early galaxies to 
argue that the dust lane in NG C 44 76 did not o riginate from the stellar component but was 
instead externally accreted. iTomita et al.l (120001 ) confirmed the existence of the regular dust 
lane using HST E(V-I) [F555W - F814W\ extinction images and showed that NGC 4476 has 
H II regions. They also noted a blue unresolved nucleus that th e y sug gested may come from 
a young star cluster(s) associated with the dust system. lYoungl (120021 ) detected CO gas in a 
relaxed disk that is co-spatial with the dust lane and has a diameter of about 27". 

Figure [7] shows 350 /xm-continuum contours of dust emission from NGC 4476 that gen- 
erally follow the direction of the ellipticity and major axis position angle of those of the CO 
disk and SDSS-extinction maps except at parts of the 350 /im-continuum at 30% flux levels. 
In the eastern part of the galaxy, the 70%, 50%, and 30% contour levels of the SDSS i-band 
and SHARC II 350 /mi (respectively of the stellar and dust components) intensities co-incide 
quite well. The SDSS extinction image in this region all extends to the 30% SHARC II 
350 yum contour level. However, in the western part of the galaxy, the correspondence of the 
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contours levels is not so good, suggesting that the dust and stellar disks are not in the same 
plane. Indeed, the SDSS-color map shows less extinction in the western part of the galaxy, 
as may be the case if the dust disk is towards the background or inner parts of the galaxy in 
that region. The 350 /zm-continuum 50% contour levels along the major axis of the galaxy 
are more extended than the CO and stellar intensities of the same relative levels. These 
results suggest that the dus t component is separate from the stellar one, as concluded by 



van Dokkum fc Franxl (119951 ) from their HST data of the central region in this galaxy. The 
350 /im-continuum 30% contour levels (which are not detected at a high signal-to-noise) are 
also more extended than the CO and extinction ones and have a spur in the southwest. This 
emission spur may be part of the dust disk in the eastern background of the galaxy. 

The H II regions in NGC 4476 (and NGC 5666) most probably obtain their fuel from the 
dust seen both in the optical extinction and submm and in turn heat the dust to produce 
the FIR-submm emission seen in these galaxies. Futher discussion of star-formation, the 
heating mechanisms, as well as the dust sources are discussed in Section [5] below and will 
subsequently be analysed in more detail in a follow-up paper that uses Spitzer mid-IR 
spectroscopy of these galaxies (Leeuw et al. in prep). 



4. The Mid-IR to Submm Spectral Energy Distribution and Luminosities 

The mid-infrared (mid-IR or MIR) flux in ellipticals comes primarily from "warm" 
galactic dust, nebulae, or the envelopes of evolved stars, while the FIR to submm flux 
comes mainly from cold dust with temperatures, T < 100 K; and, in the case of radio-bright 
sources, it can also originate from the high-frequency-radio components. On the other hand, 
flux in the optical/near-IR comes primarily from the stellar components of T > 1000 K, with 
negligible contribution to the FIR to submm spectrum. For the radio-faint sources observed 
in this program, it is assumed that the mid-IR to submm SED emission comes from a cold 
to cool (5 K< T < 100 K) thermal component plus a mid-IR power-law component that 
includes PAHs, warm (T > 100 K) dust, and stellar galactic components. 

Motivated by the general appearance of the SEDs of the elliptical galaxies and the 
available data, the SED is modeled with a power law in the mid-IR (8 to 40 /im) and a 
graybody in the far-IR (40 to 2000 //m). The Mid-IR power- law fit is not a physical model, 
but rather a mathematical tool to aid in the calculation of total mid-IR luminosity for 
the galaxies. However, the graybody fit is an attempt at physically modeling the average 
temperature, emissivity, and optical depth of the cold dust emission for the galaxies. 

The model fits are plotted in Figure [HI and their derived parameters are labeled on the 
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plots as well as listed in Table HI The power-law model is generated with a truncated Fv 
temperature distribution, and thus the shape of the mid-IR components shown in the SED 
plots. For sources which have less than two measurements in the mid-IR, the median power 
law Fv 4 is assumed. The power law for NGC4476, which has no mid-IR detection, shows 
a 3 cr upper limit in Figure Reasonable temperature and beta are derived for the graybody 
part of the spectrum for all sources except NGC 83. For that source, it is assumed that the 
dust emissivity index beta is equal to the median value of 1.8. 

Cold dust temperatures and emissivity indices are respectively determined to be between 
22 K < T < 31 K and 1.6 < (3 < 2.2 for the sample. These parameters are typical for single- 



Leeuw et al. 


2004; 


Temi et al. 


2004) 



(e.g., 

than dust temperatures of ellipticals derived from 60 and 100 fim flux ratios alone and a 
fixed (3 = 1.0 (e.g. Wiklind et al. 1995). On the other hand, the dust temperatures are also 
slightly warmer and the emissivi ty indices similar t o those of dust in the Milky Way Galaxy 



and nearby spiral galaxies (e.g.. lDunne et al.ll2000l ). 



The cold dust continuum emission presented here extends over 10 kpc for some galaxies 
in the sample, is detected at high signal to noise with SHARC II, and is generally co-spatial 
with previously imaged CO emision (see, e.g., Section [3]). Above the 50% contour levels, the 
submm dust continuum and optical stellar intensities of similar spatial resolution have the 
same general distribution, suggesting that the dust, CO, and stellar components in the central 
regions are all co-spatial. The early-merger elliptical NGC 3656, with its disturbed optical 
morphology, is an exception to this coupling of submm and optical central flux distributions. 

The current maps are not detected at high enough significance that proper dust distri- 
butions can be discerned below the 30% contour levels. If there exists any dust component 
for the sample that is cooler or has lower-level emission that perhaps extends over the entire 
galactic stellar or HI distributions (see, e.g., Section [5]), its detection will have to await the 
next-generation submm detectors such as those that will be on ALMA, Herschel, JCMT, 
or SOFIA. Combining the current data set together with anticipated Spitzer and other up- 
coming sensitive FIR-to-submm observations and analysing the SED of mapping data from 
various FIR-to-submm bands should put tighter constraints on the physical parameters of 
cold dust in the sample. 

MIR-to-submm luminosities (£o.o4-imm an d is/im-imm) that were obtained by integrat- 
ing below the model fits shown in Figure [H] are listed in Table and these can be compared 
to the optical luminosities (£b) that are also listed in Table The percentage excess MIR- 
to-submm luminosities over the total luminosities (LB-imm) are determined and also listed 
in Table [5j For most of the galaxies, the MIR-to-submm excess luminosities are about 20% 
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to 30% of the total luminosities, with about 70% or more of the excess emitted at FIR to 
mm wavebands as thermal, cold dust emission (see Figure [8] and Table [5]). The respective 
early-age and intermediate-age galaxies NGC 3656 and NGC 5666 have higher than average 
MIR-to-submm excesses of the total luminosities ( > 50%) for the sample (see Table [5]), 
suggestive of more or denser dust and star-formation or, at least, fairly high MIR-to-submm 
emission in comparison to total emission in these more recent merger- remnant/dust-accretion 
galaxies. 



5. Estimates of Star Formation Rates and Other Luminosity Sources 

The peak 350 /im optical depths given in Table [6] were determined using the 350 /im 
continuum fluxes and intrinsic 350 /x m sizes respectively presented in Tables |2] and [3] for each 



galaxy. According to iDraind (120031 ) . the corresponding peak optical extinctions are about 
14,000 times larger than the sub-mm optical depth values, implying peak visual extinctions 
{Ay) greater than 2 mag for all the galaxies in the sample. Therefore, the dust measured at 
350 yum for these galaxies is dense enough to absorb almost all the UV or optical radiation 
emitted by any massive star embedded within it. 

If the far-IR luminosity, L(FIR), for each galaxy in our sample is due to the reprocessing 
of radiation from embedded, newly formed stars, then an estimate can be made of the current 
Star Formation Rate, SFR, in each galaxy using the relation given in Kennicutt (1998) 

SFR(M Q /yr) = 1.75 x 10" 10 L Q (FIR), (1) 

where L(FIR) includes radiation from 8 through to 1000 /im. Table [5] lists the estimate of 
L(FIR) for each galaxy based on the fits to the data shown in Figure[Hl Except for NGC 4476, 
Equation 1 yields SFR estimates ranging from 1 to 3 M Q /yr for the galaxies in the sample 
(see Table [6]) . These rates are similar to the SFR in the Milky Way, a spiral galaxy, but are 
about 10 times those for 'normal' elliptical galaxies. NGC 4476 has a SFR estimate more 
in line with a normal elliptical galaxy. However, NGC 4476 is exceptional in this sample in 
that it is a dwarf elliptical near the center of the Virgo cluster. 

For the most part, the FIR luminosity given in Table [5] for each galaxy is composed of 
power radiating at wavelengths longer than 40 /im. Based on the SED's shown in Figure [HJ 
the majority of this FIR power, therefore, could be coincident with the 350 /im emission, 
which seems to be embedded in rotating CO "disks" . Thus, the FIR data and our 350 /im 
data imply for most of the galaxies in our sample, the presence of mini-starbursts inside 
central rotating disks of gas and dust. 



There is another more indirect way to estimate SFR, by using the empirical relationship 
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from iKennicuttl (119981 ) between SFR surface density, SD/^fk), an d the molecular gas surface 



density, ST> (Gas) : 



,SD (SFR) (M /yr/kpc 2 ) = 2.5 x 10-\SD {Gas) (M Q /pc 2 )) 1A . (2) 



Estimates of the SFR surface densities at the peaks of the gas distributions in these galaxies 
have been derived based on the total molecular gas masses (Table CP and the intrinsic sizes 
of the CO distributions in these galaxies (Table [3]). The product of the SFR surface densities 
with the intrinsic areas of the CO distributions yield the total SFRs for the galaxies and 
these are listed in Table [61 

A comparison between the SFR values estimated using the two methods given above 
would imply that NGC 759 and, especially, NGC 3656 have the potential (in gas) for consid- 
erably more star formation than is presently being observed in the FIR, but that the other 
galaxies in our sample are using their gas reservoirs at about the "normal" rate to form stars. 

The discussion above assumes the FIR ex cess from these galaxies orig inates from star 



formation; however, as noted in the Section 1, iMathews fc Brighentil (120031 ) argue that the 
FIR excess could be due to the thermal energy released by the hot ISM as it cools through 
thermal electron collisions with dust grains. In this scenario, the central cold dust seen in 
these elliptical galaxies originates from the settling of the cooled hot dusty gas from mass 
losing red giants distributed throughout each of the elliptical galaxies, rather than fr o m the 



gas and dust from recent galaxy merger or accretion events. IMathews &: Brighentil (120031 ) 
estimated the FIR luminosity radiated from the dust grains involved in the cooling of the hot 
gas of about 10 7 K for a given total stellar mass-loss rate, M/M yr _1 , to be approximately 

L (FIR) (L ) ~ 5.1 x lO^T/lO^XM/Moyr- 1 ). (3) 

Substituting the FIR luminosity listed for each galaxy in Table|5]gives mass-loss rates ranging 
from 15 to 340 M Q yr _1 for this set of galaxies. These rates are much higher than those for 
giant elliptical galaxies, which are normally ~ lM yr _1 . The larger mass-loss rates seem 
consistent with the larger than normal dust masses at the center of these galaxies as implied 
by our 350 yum data; however, the optical/near-IR luminosity for each galaxy in the sample 
presented here is similar to normal giant elliptical galaxy optical/near-IR luminosities, so 
normal stellar mass-loss rates would be expected for these galaxies. Hence, although both 
the FIR luminosity and the 350 fim dust mass estimates are consistent with large stellar 
mass- loss rates for the galaxies in this sample, the optical/near-IR luminosities for these 
galaxies are not consistent with such a stellar mass-loss increase. 



In summary, the calculations given in this section favor the FIR luminosity of these 
galaxies resulting primarily from cold dust in the central currently star formation regions 
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that were probably induced and supported by accretion or merger events involving a gas- 
or dust-rich galaxy or satelite, rather than coming from hot ISM cooling in collisions with 
dusty gas from stellar mass-loss. However, a small FIR luminosity contribution from the hot 
ISM cooling, the stellar, or any other dilute galactic radiation field is not excluded. 



Dust Masses and Gas-to-Dust Mass Ratios 



Following iHildebrandl (119831 ). the mass of dust, Md, that emits mid-IR to submm con- 



tinuum can be estimated from 

F D 2 

Md = I Wf TV & 

where F v is the measured flux density at frequency u, D is the distance to the source, 
B(y,T) the Planck function and k d = 3Q u /4ap the grain mass absorption coefficient where 
a and p are respecti yely the grain radius and density. A recently updated value of ^ 50 ^ m = 



1.92 x 10 1 m 2 kg 1 (lDrainell2003l ). as for Galactic dust (cf. above), is assumed, yielding dust 
masses that range from ~9x 10 5 to ~ 2.5 x 10 7 M Q for dust with temperatures in the range 
22 K< T < 30 K (see Figure E and Table 0}. 

NGC4476 has the lowest dust mass in the sample. The relative p aucity of dust i n this 



galaxy is not surprising as NGC4476 is a dwarf elliptical. However, iLucero et al.l (120051 ) 
observed that this galaxy, which lies within the Virgo Cluster, has been stripped of its HI, 
maybe through ram-pressure; therefore, the presence of any compact distribution of CO and 
cold dust at all implies that the H I is more vulnerable to ram-pressure stripping than is the 
denser molecular gas and the dust associated with it within this cluster galaxy. 

Table H] shows that the sample has diverse H2-mass to dust-mass ratios that are mostly 
between 120 < M(H 2 )/Md ust < 160, and in the quiescent galaxy NGC807 and early-merger 
remnant NGC 3656 respectively range from ~ 60 (about half the average) to ~ 310 (about 
twice the average). The ratios suggest that, in the central galactic regions where the H 2 and 
dust masses are respectively measured from the resolved CO and our CSO continuum mea- 
surements, some of the sample galaxies are relatively under-massive in H 2 and over-massive 
in dust content, while the others are over- massive in H 2 and under- massive in dust content. 
The implications of these mass ratios and their possible association with star-formation or, 
at least, gas-phases in the central regions of these galaxies will be investigated further in 
our upcoming paper that uses Spitzer IRS data to study the MIR spectral evolution and/or 
star- formation properties of these galaxies. 

The total cold-gas-mass (i.e., (H 2 + HI)-mass) to dust-mass ratios for the sample are, 
however, less diverse and in a narrower range of 230 < M(H 2 + HI)/Md us t < 400, with 
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NGC807 and NGC3656 at much closer ratios of 300 and 387, respectively (see Table Hj). At 
face value, it may appear that the total cold-gas-mass to dust-mass ratio for the sample is 



about 3 times the value for the Milky Way Galaxy (i.e. , M(H? + HP/M , 



Lake et al. 1987: 



"St, 



100). However , 



Balcells et al.l l200lh . 



in sample objects where HI has been resolved (e.g., 
the atomic gas has been observed to be much more extended than the inner galactic regions 
where our observations have mapped dust emission and shown it to be associated with 
molecular gas. A cold and extended dust component may yet be detected that is associated 
with the atomic hydrogen, and perhaps massive enough to bring the total cold gas-to-dust 
mass ratio for the sample in line with the Milky Way value. 



7. Summary of Results and Conclusions 

SHARCII continuum observat ions at 350 /im have been presented o f seven elliptical 



galaxies with known CO-gas disks (IWiklind et al.l If 9971 ; lYoung 



2002 . 120051 ) . The SHARC II 



mapping provides the first clearly resolved FIR-submm continuum emission from cold dust 
(with temperatures ~ 31 K > T > 23 K) of any elliptical galaxy > 40Mpc. Calculations of 
the measured FIR excess show that the dust's most likely and dominant heating source is 
not dilute stellar radiation or cooling flows but rather star-formation, that could have been 
triggered by an accretion or merger event and fueled by dust-rich material that has settled 
in a dense region co-spatial with the central CO-gas disks. The dust is detected even in two 
cluster ellipticals that are deficient in HI, showing that, unlike the HI, cold dust and CO in 
ellipticals can survive among hot X-ray gas, even in galaxy clusters. 

Above the 50% contour levels, where the submm detections are of high significance, 
the submm dust continuum and optical stellar intensities of similar spatial resolution have 
the same general distribution, suggesting that the dust, CO, and stellar components in the 
central regions are co-spatial. The early-merger elliptical NGC3656, with its dense optical 
dust lane that is thought to be in a highly inclined warped disk, is an exception to this 
coupling of submm and optical central flux distributions. Below the 50% contour levels, and 
especially below the 30% levels, the current maps are not detected at high enough significance 
that proper dust distributions can be discerned. 

No dust cooler than 20 K, distributed outside the CO disks, or co-spatial with and 
heated by the entire, dilute stellar optical galaxy (or very extended HI) is currently evident. 
These observations clearly show that the FIR-submm emission in our sample ellipticals is 
primarily from cold dust that is associated with star-formation located in the central parts 
of these galaxies. The results support the hypothesis that if ellipticals are formed through 
mergers of gas- rich spirals, their starbursts (1) are necessary to produce the high phase- 
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space density seen in the center of ellipticals and (2) would power the FIR-bright stage as 
the merger relaxes. Combining the current data set together with anticipated Spitzer and 
other upcoming sensitive FIR-to-submm mapping observations may detect or rule-out the 
existence of larger-scale, diffusely distributed cold dust with an average temperature that is 
perhaps < 20 K (and maybe associated with HI). 

Assuming the grain mass a bsorption coefficient /^ 50 ^ m = 1.92 x 10 -1 m 2 kg -1 , as for 



Galactic-like dust (iDraind 120031 ) . the measured submm fluxes and dust temperatures yield 
dust masses for the sample that range from ~ 9 x 1O 5 M to ~ 2 x 10 7 M Q . The sample has 
diverse H 2 -mass to dust-mass ratios that cover 60 < M(H 2 )/Md us t < 310, suggesting that in 
the central galactic regions, where the H2 and dust masses are respectively measured from the 
resolved CO and our CSO/SHARCII continuum measurements, some of the sample galaxies 
are relatively under-massive in H 2 and over-massive in dust content, while others are the 
opposite. The possible implications of these mass ratios on star- formation or, at least, gas- 
phases in the central regions of these galaxies will be investigated further in our upcoming 
paper that uses Spitzer IRS data to study the detailed MIR spectral and/or star-formation 
properties of these galaxies. 

The total cold-gas-mass to dust-mass ratios for the sample are in a narrower range of 
230 < M(H 2 + HI)/Mdust ;$ 400, with a total cold-gas-mass to dust-mass ratio for the 
sample that, at face value, may appear to be about 3 times the value for the Milky Way 
Galaxy (i.e., ~ 100)- However, in sa mple objects where atomic gas has been resolved (e.g., 



Lake et al.lll987l ; iBalcells et al.ll200ll ). the HI has been observed to be much more extended 
than the inner galactic regions where our observations have mapped dust emission that is 
associated with molecular gas. A cold and extended dust component may yet be detected 
that is associated with the atomic hydrogen, and perhaps massive enough to bring the total 
cold gas-to-dust mass ratio for the sample in line with the Milky Way value. As noted above, 
this may be oberved with future more sensitive FIR/submm detectors such as those that 
will be on SOFIA, Herschel, JCMT, or ALMA. 



Part of this work was done while Lerothodi L. Leeuw (LLL) was at the University of 
Chicago and partially supported by a Chicago 2004 NASA Mini- Award Grant and NSF grant 
AST-0505124 and at NASA Ames Research Center. LLL finalized this research and wrote 
this paper while supported by the South Africa Square Kilometer Array Postdoctoral Bursary 
at Rhodes University. Research at the Caltech Submillimeter Observatory is supported by 
NSF grant AST-0540882. 

This publication makes use of data products from the second Digital Sky Survey (DSS2), 
that was based on photographic data obtained using the Oschin Schmidt Telescope on Palo- 



-17- 



mar Mountain. The Palomar Observatory Sky Survey was funded by the National Geo- 
graphic Society. The Oschin Schmidt Telescope is operated by the California Institute of 
Technology and Palomar Observatory. The plates were processed into the present com- 
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Table 1: Publi s hed optical and molecular-gas propert i es of the ell i ptica l s (ad opted from 



Murray et al.1 l200d : iTomita et al.1 12000 : IWiklind et al.1 11995L Il997l : lYound 12002 . 12003 . and refer- 
ences therein) 



Source 


Type 


Environment, Notables, & 


M(H 2 ) 


M(HI) 


D 


L B 


(B~ 






Merger or Accretion Stage 


[10 8 


[10 8 


[M 


[10 9 


V)e 








M ] 


Mo] 


pc] 


L©] 


mag 


NGC 3656 


Ep 


Early-Age Major Merger/Acer.? 


47 


11 


45 


16 




NGC 5666 


cE2/SO 


Intermediate Merger / Acer. ? 


5.7 


10 


35 


6 


0.86 


NGC 83 


E0 


Group Member, Late Merger/Acer.? 


20 




85 


44 


1.12 


NGC 759 


E0/1 


Cluster A262, Late Merger /Acer.? 


24 


<24 


66 


34 


1.05 


NGC 807 


E3 


Field, Quiescent / Very-Late Acer.? 


14 


61 


64 


31 


0.97 


UGC 1503 


El 


Field, Quiescent / Very-Late Acer.? 


18 


16 


69 


16 




NGC 4476 


dE5p/SO 


Virgo Cluster, Dwarf Merger/Acer.? 


1.1 


<1 


18 


3 


0.85 



Table 2: SHARC II & IRAS continuum fluxes for the selected elliptical galaxies 



Source 


Peak @350/mi 


Total 


Total 


Total 


Total 


Total 




(mJy/9" beam) 


350/iin (Jy) 


100//m (Jy) 


60/im (Jy) 


25/im (Jy) 


12/im(Jy) 


NGC 3656 


368 ± 73 


0.69 ±0.14 


6.58 ±0.68 


2.54 ±0.13 


0.32 ±0.03 


0.14 ±0.03 


NGC 5666 


150 ±30 


0.52 ±0.10 


3.98 ±0.40 


1.99 ±0.10 


0.16 ±0.04 


0.12 ±0.04 


NGC 83 


191 ±38 


0.25 ±0.05 


2.15 ±0.28 


0.34 ±0.09 


< 0.03 


0.06 ±0.02 


NGC 759 


71 ± 14 


0.42 ±0.08 


2.36 ±0.28 


0.85 ±0.06 


0.07 ±0.02 


0.07 ±0.03 


NGC 807 


63 ± 12 


0.45 ±0.09 


1.83 ±0.22 


0.41 ±0.03 


0.12 ±0.02 


0.12 ±0.03 


UGC 1503 


57 ±11 


0.24 ±0.05 


1.43 ±0.20 


0.40 ±0.03 


0.07 ±0.03 


< 0.03 


NGC 4476 


117 ±23 


0.28 ±0.06 


1.84 ±0.21 


0.66 ±0.05 


< 0.04 


< 0.05 
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Table 3: Intrinsic Galactic sizes of the stellar, molecular gas, and dust distributions at the 50% 
contour levels of the respective peak surface brightness (The listed SDSS i-band size for NGC 3656 
is at the 70% contour level. See text for further details.) 

Source @50% 350 /xm Flux @50% CO Flux @50% SDSS i-Band Flux 

or * @50% DSS2 Red Flux 





arcs x arcs 


kpc x kpc 


arcs x arcs 


kpc x kpc 


arcs x arcs 


kpc x kpc 


NGC 3656 


6.5 x 12.7 


1.5 x 2.8 


2.6 x 13.2 


0.5 x 2.8 


11.2 x 17.3 


2.4x 3.8 


NGC 5666 


10.3 x 11.2 


1.8 x 2.0 


8.2 x 11.1 


1.3 x 1.8 


12.5 x 15 


2.1 x 2.5 


NGC 83 


4.1 x 5.0 


2.0 x 1.8 


3.8 x 5.2 


1.5 x 2.0 


saturated 


saturated 


NGC 759 


5.0 x 5.0 


1.5 x 1.5 


3.7 x 4.5 


1.0 x 1.5 


saturated 


saturated 


NGC 807 


13.7 x 30.4 


4.3 x 9.3 


11.7 x 20.3 


3.5 x 6.3 


*20 x 32 


*6.3 x 9.8 


UGC 1503 


12.1 x 21.4 


4.0 x 7.0 


9.4 x 15.2 


3.0 x 5.0 


*15.3 x 18.9 


*5.0 x 6.3 


NGC 4476 


4.1 x 19.1 


0.3 x 1.5 


3.0 x 13.4 


0.3 x 1.0 


2.0 x 9.8 


0.2 x 0.9 



Table 4: Mid-IR to SHARC II-derived cold-dust parameters and related gas properties 



Source 


T 
[K] 


P 


M(d) 
[1O 6 M ] 


M(U 2 ) 
M(d) 


M(d) 


NGC 3656 


25.4 


2.2 


15 


313 


387 


NGC 5666 


30.7 


1.6 


4.8 


119 


327 


NGC 83 


25.8 


(1.8) 


19 


105 




NGC 759 


29.1 


1.5 


15 


160 


<320 


NGC 807 


22.6 


1.9 


25 


56 


300 


UGC 1503 


25.1 


1.9 


12 


150 


283 


NGC 4476 


27.2 


1.8 


0.9 


122 


<233 



Table 5: Mid-IR to mm luminosity excess 



Source 


r 0.04mm 
to 1mm 


r 8(im 
to 1mm 


T B to 
1mm 


0/ J B to 
/0L lmin 


(V T B to 
/0 -klmm 




[10 9 


[10 9 


[10 9 


as 


as 




L ] 


L ] 


L ] 


-^to 1mm 


j 0.04mm 
to 1mm 


NGC 3656 


14.0 


17.8 


33.8 


53 


41 


NGC 5666 


5.4 


6.8 


12.8 


53 


42 


NGC 83 


11.7 


13.8 


57.8 


24 


20 


NGC 759 


10.9 


13.8 


47.8 


29 


23 


NGC 807 


8.1 


12.6 


43.6 


29 


18 


UGC 1503 


6.9 


8.3 


24.3 


34 


28 


NGC 4476 


0.6 


0.7 


3.7 


19 


16 
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Table 6: The optical depth and total FIR and CO star-formation rate estimates 



Source 


350 fim 


A v 


Total 


Total 




Peak 




FIR SFR 


CO SFR 




^optical 




M /yr 


M /yr 


NGC3656 


0.00119 


17 


3.04 


26.89 


NGC5666 


0.00050 


7 


1.17 


1.17 


NGC83 


0.00171 


24 


2.65 


6.13 


NGC 759 


0.00206 


29 


2.41 


10.38 


NGC807 


0.00018 


25 


2.17 


1.73 


UGC1503 


0.00015 


2 


1.44 


2.83 


NGC 4476 


0.00047 


7 


0.14 


0.28 
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Fig. 1.— Top Left: White contours of the total integrated CO (1-0) intensity map of 
NGC3656, with a measured 10%-contour-level diameter of 34", are overlayed on the grey 
scale and black contou r s from the red part of the optical Digitized Sky Survey 2 (DSS2), 
as presented by lYoungl (120021 ). The white contours are in units of -5%, -2%, 2%, 5 %, 10%, 



20%, 30%, 50%, 70%, and 90% of 81.1Jy beam^kms- 1 = 4.7 x 10 22 cm" 2 CO peak ((Young 



20021 ). Top Right: White contours of the CSO/SHARCII 350 /im continuum smoothed to 
10" are overlayed on intensity-grey scale and black contours of the Sloan Digital Sky Survey 
(SDSS) i-band image of NGC3656 smoothed to 10". The SDSS i-band and submm contours 
are respectively ~ 57%, 58%, 59%, 60%, 70%, and 90% and 10%, 30%, 50%, 70%, and 90% 
of the maximum pixel values on the maps and are intended to display areas that may be co- 
spatial with the CO as shown in the adjacent figure. Bottom left: A zoom- in of the top-right 
figure. Bottom right: White contours of the CSO/SHARCII 350 continuum smoothed 
to ~ 10" are overlayed on grey scale and black contours of the SDSS i- minus g-band images 
smoothed to 2". The SDSS contours are on a linear scale, and the submm contours are as 
above. In the SDSS i- minus g-band maps, darker pixels have redder colors. 
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Fig. 2.— Top Left: White contours of the total integrated CO (1-0) intensity map of 
NGC5666, with a measured 10%-contour-level diameter of 26", are overlayed on the grey 
scale and black contour s from the red part of the optical Digitized Sky Survey 2 (DSS2), 
as presented by lYoungl (120021 ) . The white contours are in units of -10%, -5%, 5% , 10%, 
20% 30%, 50%, 70%, and 90% of 21.3 Jy beanrtms- 1 = 7.5 x 10 21 cm~ 2 CO peak jYoung 
20021 ). Top Right: White contours of the CSO/SHARCII 350 /im continuum smoothed to 
10" are overlayed on intensity-grey scale and black contours of the Sloan Digital Sky Survey 
(SDSS) i-band image of NGC5666 smoothed to 10". The SDSS i-band and submm contours 
are respectively ~ 40%, 41%, 42%, 45%, 50%, 70%, and 90% and 20%, 30%, 50%, 70%, 
and 90% of the maximum pixel values on the maps and are intended to display areas that 
may be co-spatial with the CO as shown in the adjacent figure. Bottom left: A zoom- in of 
the top-right plot. Bottom right: White contours of the CSO/SHARCII 350 /im continuum 
smoothed to ~ 10" are overlayed on grey scale and black contours of the SDSS i- minus 
g-band images smoothed to 2". The SDSS contours are on a linear scale, and the submm 
contours are as above. In the SDSS i- minus g-band maps, darker pixels have redder colors. 
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Fig. 3. — Left: White contours of the total integrated CO(l-O) intensity map of NGC83, 
with a measured 10%-contour-level diameter of 10", are overlayed on the grey scale and 
bla ck cont o urs fr om the red part of the optical Digitized Sky Survey 2 (DSS2), as presented 
bv lYoungi (boosl ). The white contours are in units of -20%, -10%, 1 0%, 20%, 30 %, 50%, 
70%, and 90% of the 12.7 Jy beam^kms- 1 = 8.8 x 10 21 cm~ 2 CO peak jYoiinjboOfih . Right- 
White contours of the CSO/SHARCII 350 /im continuum smoothed to 10" are overlayed on 
intensity-grey scale and black contours of the DSS2 blue image of NGC 83. The DSS2-blue 
and submm contours are respectively 50%, 70%, and 90% and 20%, 30%, 50%, 70%, and 
90% of the maximum pixel values on the maps and are intended to display areas that may 
be co-spatial with the CO as shown in the adjacent figure. The DSS2 image of NGC 83 is 
saturated and plotted only as a guide of the possible optical extent of this galaxy. 



- 25 - 



36 20 40 - 



36 20 30 - 



36 20 25 



36 20 20 h_ 





£ 36:20:30 
Q 



I 1 1 1 I 1 1 1 I 



I 1 1 1 I 1 1 1 I 1 '.' I 1 1 1 I 1 1 1 I 1 




51.2 50.8 50.4 1:57:50.0 49.6 
Right ascension 



Fig. 4.— Top Left: Total integrated CO(l-O) i ntensity map of NG C 759 (grey scale and 
black contours) that is adopted from the paper by lWiklind et al.l (119971 ). The diameter of the 
CO, black contou r level spacing, and rm s noise in the map are respectively 5", 0.8 Jy kms -1 , 



and 0.8 Jy kms" 1 IWiklind et all (119971 ) Top Right: White contours of the CSO/SHARCII 



350 yum continuum smoothed to 10" are overlayed on intensity-grey scale and black contours 
of the Digitized Sky Survey 2 (DSS2) blue image of NGC 759. The DSS2-blue and submm 
contours are respectively 30%, 50%, 70%, and 90% and 20%, 30%, 50%, 70%, and 90% of the 
maximum pixel values on the maps and are intended to display areas that may be co-spatial 
with the CO as shown in the adjacent figure. 
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Fig. 5.— Top Left: White contours of the total integrated CO (1-0) intensity map of 
UGC 1503, with a measured 10%-contour-level diameter of 30", are overlayed on the grey 
scale and black contou r s from the red part of the optical Digitized Sky Survey 2 (DSS2), 
as presented by lYoungi (120021 ) . The white contours are in units of -20%, -10%, 10%, 20%, 
30%, 50%, 70%, an d 90% of 6.3 Jy beam^kms- 1 = 3.9 x 10 21 cm- 2 CO integrated intensity 
peak (lYoungi |2002| ). Top Right: White contours of the CSO/SHARCII 350 /im continuum 
smoothed to 10" are overlayed on intensity-grey scale and black contours of the DSS2 blue 
image of UGC 1503. The DSS2-blue and submm contours are respectively 20%, 30%, 50%, 
70%, and 90% and 20%, 30%, 50%, 70%, and 90% of the maximum pixel values on the 
maps and are intended to display areas that may be co-spatial with the CO as shown in the 
adjacent figure. 
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Fig. 6.— Top Left: White contours of the total integrated CO (1-0) intensity map of 
NGC807, with a measured 10%-contour-level diameter of 40", are overlayed on the grey 
scale and black contou r s from the red part of the optical Digitized Sky Survey 2 (DSS2), 
as presented by lYoung) (120021 ) . The white contours are in units of -20%, -10 %, 10%, 20% , 
30%, 50%, 70%, and 90% of 7.6 Jy beam^kms- 1 = 2.6 x 10 22 cm~ 2 CO peak jYounejl2002h . 
Top Right: White contours of the CSO/SHARCII 350 /im continuum smoothed to 10" are 
overlayed on intensity-grey scale and black contours of the DSS2 blue image of NGC807. 
The DSS2-blue and submm contours are respectively 20%, 30%, 50%, 70%, and 90% and 
20%, 30%, 50%, 70%, and 90% of the maximum pixel values on the maps and are intended 
to display areas that may be co-spatial with the CO as shown in the adjacent figure. 
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Fig. 7.— Top Left: White contours of the total integrated CO (1-0) intensity map of 
NGC4476, with a measured 10%-contour-level diameter of 27", are overlayed on the grey 
scale and black contou r s from the red part of the optical Digitized Sky Survey 2 (DSS2), 
as presented by lYoungi (120021 ) . The white contours are in units of -20%, -10 %, 10%, 20% , 
30%, 50%, 70%, and 90% of 12.4 Jy beam^kms- 1 = 7.3 x 10 21 cm" 2 CO peak jYoungll2002h . 
Top Right: White contours of the CSO/SHARCII 350 /im continuum smoothed to 10" are 
overlayed on intensity-grey scale and black contours of the Sloan Digital Sky Survey (SDSS) 
i-band image of NGC4476 smoothed to 10". The SDSS i-band and submm contours are 
respectively ~ 31%, 32%, 33%, 35%, 40%, 50%, 70%, and 90% and 30%, 50%, 70%, and 
90% of the maximum pixel values on the maps and are intended to display areas that may 
be co-spatial with the CO as shown in the adjacent figure. Bottom left: A zoom in of the 
top-rright figure. Bottom right: White contours of the CSO/SHARCII 350 /im continuum 
smoothed to ~ 10" are overlayed on grey scale and black contours of the SDSS i- minus 
g-band images smoothed to 2". The SDSS contours are on a linear scale, and the submm 
contours are as above. In the SDSS i- minus g-band maps, darker pixels have redder colors. 
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Fig. 8. — The CSO/SHARCII 350 yum continuum integrated fluxes and mid-IR to submm 
data obtained from public archives and literature (see Table [2]) for the sample galaxies are 
plotted and fitted with a composite model of a single-temperature greybody plus mid-IR 
power-law (see Section Hj). 



